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We develop the synchronization theory of microwave induced zero-resistance states (ZRS) for 
two-dimensional electron gas in a magnetic field. In this theory the dissipative effects lead to syn- 
chronization of cyclotron phase with driving microwave phase at certain resonant ratios between 
microwave and cyclotron frequencies. This synchronization produces stabilization of electron trans- 
port along edge channels and at the same time it gives suppression of dissipative scattering on local 
impurities and dissipative conductivity in the bulk, thus creating the ZRS phases at that frequency 
ratios. The electron dynamics along edge and around circular disk impurity is well described by 
the Chirikov standard map. The theoretical analysis is based on extensive numerical simulations of 
classical electron transport in a strongly nonlinear regime. We also discuss the value of activation en- 
ergy obtained in our model and the experimental signatures that could establish the synchronization 
origin of ZRS. 



PACS numbers: 73.40.-c,05.45.-a,72.20.My 

I. INTRODUCTION 

The experiments on resistivity of high mobility two- 
dimensional electron gas (2DEG) in presence of a rela- 
tively weak magnetic field and microwave radiation led to 
a discovery of striking Zero-Resistance States (ZRS) in- 
duced by a microwave field by Mani et al. and Zudov 
et al. [2]. Other experimental groups also found the mi- 
crowave induced ZRS in various 2DEG samples (see e.g. 
[3H5]). A similar behavior of resistivity is also observed 
for electrons on a surface of liquid helium in presence 
of magnetic and microwave fields [5J [7]. These experi- 
mental results obtained with different systems stress the 
generic nature of ZRS. Various theoretical explications 
for this striking phenomenon have been proposed during 
the decade after the first experiments [J [5] . An overview 
of experimental and theoretical results is give in the re- 
cent review ^ . 

In our opinion the most intriguing feature of ZRS is 
their almost periodic structure as a function of the ra- 
tio j = 10 1 10 c between the microwave frequency lo and 
cyclotron frequency lo^ = eB/mc. (in the following we 
are using units with electron charge e and mass m equal 
to unity). Indeed, a Hamiltonian of electron in a mag- 
netic field is equivalent to an oscillator, it has a magneto- 
plasmon resonance at j = 1 but in a linear oscillator there 
are no matrix elements at j — 2, 3, ... and hence a rela- 
tively weak microwave field is not expected to affect elec- 
tron dynamics and resistivity properties of transport. Of 
course, one can argue that impurities can generate har- 
monics being resonant at high j > 1 but ZRS is observed 
only in high mobility samples and thus the density of im- 
purities is expected to be rather low. It is also important 
to note that ZRS appears at high Landau levels z/ ~ 50 
so that a semiclassical analysis of the phenomenon seems 
to be rather relevant. 

In this work we develop the theoretical approach pro- 



posed in [S]. This approach argues that impurities pro- 
duce only smooth potential variations inside a bulk of 
a sample so that ZRS at high j appear from the orbits 
moving along sharp sample boundaries. It is shown [5] 
that collisions with boundaries naturally generate high 
harmonics and that a moderate microwave field gives 
stabilization of edge channel transport of electrons in a 
vicinity of j « = 1 1/4, 2 + 1/4, 3 1/4... pro- 
ducing at these j a resistance going to zero with increas- 
ing microwave power. This theory is based on classical 
dynamics of electrons along a sharp edge. The treat- 
ment of relaxation processes is modeled in a phenomeno- 
logical way by a dissipative term in the Newton equa- 
tions. Additional noise term in the dynamical equations 
takes into account thermal fluctuations. The dissipation 
leads to synchronization of cyclotron phase with a phase 
of microwave field producing stabilization of edge trans- 
port along the edges in a vicinity of resonant jV values. 
Thus, according to the edge stabilization theory [9 the 
ZRS phase is related to a universal synchronization phe- 
nomenon which is a well established concept in nonlinear 
sciences [TU] . 

While the description of edge transport stabilization 
[9] captures a number of important features observed in 
ZRS experiments it assumes that the contribution of bulk 
orbits in transport is negligibly small. This assumption 
is justified for smooth potential variations inside the bulk 
of a sample. However, a presence of isolated small scale 
scatterers inside the bulk combined with a smooth po- 
tential component can significantly affect the transport 
properties of electrons (see e.g. PT). Also the majority 
of theoretical explanations of ZRS phenomenon considers 
only a contribution of scattering in a bulk [5]. Thus it 
is necessary to analyze how a scattering on a single im- 
purity is affected by a combined action of magnetic and 
microwave fields. In this work we perform such an anal- 
ysis modeling impurity by a rigid circular disk of finite 
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radius. We show that the dynamics in a vicinity of disk 
has significant similarities with dynamics of orbits along 
a sharp edge leading to appearance of ZRS type features 
in a resistivity dependence on j. 

The paper is composed as follows: in Section II we 
discuss the dynamics in edge vicinity, in Section III wc 
analyze scattering on a single disk, in Section IV we study 
scattering on many disks when their density is low, here 
we determine the resistivity dependence on j and other 
system parameters, physical scales of ZRS effects are ana- 
lyzed in Section V, effects of two microwave driving fields 
and other theory predictions are considered in Section VI, 
discussion of the results is given in Section VII. 

We study various models which we list here for a reader 
convenience: wall model described by the Newton equa- 
tions ([I]), ([2]) with microwave field polarization perpen- 
dicular to the wall (model (Wl) equivalent to model 
(1) in |9]); the Chirikov standard map description ([s]) 
of the wall model dynamics called model (W2) (equiv- 
alent to model (2) in [S] at parameter p = 1); the sin- 
gle disk model with radial microwave field called model 
(DRl); the Chirikov standard map description ^ of 
model (DRl) called model (DR2) (here Vy — Vr in (j3|, 
p > 1); the model of a single disk in a linearly polar- 
ized microwave field and static electric field called model 
(Dl); the model of transport in a system with many disks 
called model (D2) which extends the model (Dl); exten- 
sion of model (D2) with disk roughness and dissipation in 
space called model (D3); the wall model (W2) extended 
to two microwave fields is called model (W3). 

II. DYNAMICS IN EDGE VICINITY 

We remind first the approach developed in [9j. Here, 
the classical electron dynamics is considered in a proxim- 
ity of the Fermi surface and in a vicinity of sample edge 
modeled as a specular wall. The motion is described by 
Newton equations 

dv/dt = X V + Loecosijjt — 7(u)v + lec + Is (1) 

where a dimensionless vector e — eE/^muvp) describes 
microwave driving field E. Here an electron veloc- 
ity V is measured in units of Fermi velocity vp and 
j{v) — 7o(|vp — 1) describes a relaxation processes to 
the Fermi surface. We also use the dimensionless am- 
plitude of velocity oscillations induced by a microwave 
field e — e\'E\/{mujvp). As in [9 , in the following we use 
units with vp = 1. The last two terms I^c and Ig in ([T]) 
account for elastic collisions with the wall and small an- 
gle scattering. Disorder scattering is modeled as random 
rotations of v by small angles in the interval ±ai with 
Poissonian distribution over time interval = l/w. The 
amplitude of noise is assumed to be relatively small so 
that the mean free path ie is much larger than the cy- 
clotron radius rc — vp/cUc- We note that the dissipative 
term is also known as a Gaussian thermostat [TJ or as a 
Landau-Stuart dissipation [lOj. The dynamical evolution 



described by Eq. ([ij is simulated numerically using the 
Runge-Kutta method. Following j9j we call this system 
model (Wl) (equivalent to model (1) in [S]). 

We note that for typical experimental ZRS parameters 
we have: electron density Ue — 3.5 • lO^^cm"^, effec- 
tive electron mass m — 0.065me, microwave frequency 
/ = a;/27r = 50GHz, Fermi energy Ep = mvp/2 = 
-ktlJt? /ra — 0.01289^, corresponding to Ep/ks = 
14:9.5Kelvin, with Fermi velocity vp — 2.641 • lO^cm/s. 
At such a frequency the cyclotron resonance w = Wc = 
eB/mc takes place at i? = O.llGlTesZa with the cy- 
clotron radius Vc = vpjiOc = 0.8873/xm. At such a 
magnetic field we have the energy spacing between Lan- 
dau levels huj = hujc — 0.2067toF = 2AQK ■ ks corre- 
sponding to a Landau level v = Ep/fiUc ~ 62. For a 
microwave field strength E = IV /cm we have the pa- 
rameter e — eE/{muJVp) — 0.003261. With these phys- 
ical values of system parameters we can always recover 
the physical quantities from our dimensionless units with 
m ^ e ~ Vp = 1. 

Examples of orbits running along the edge of specular 
wall are given in [9] (see Fig.l there). A microwave field 
creates resonances between the microwave frequency uj 
and a frequency of nonlinear oscillations of orbits col- 
liding with the wall. Due to a specular nature of this 
collisions the electron motion has high harmonics of cy- 
clotron frequency that leads to appearance of resonances 
around j — I, 2,3,4... (there is an additional shift of ap- 
proximate value 1 /4 to jr values due to a finite width of 
nonlinear resonance). 

To characterize the dynamical motion it is useful to 
construct the Poincare section following the standard 
methods of nonlinear systems [131 HI] ■ We consider the 
Hamiltonian case at 70 = in absence of noise. Also 
we choose a linear polarized microwave field being per- 
pendicular to the wall which is going along cc-axis (same 
geometry as in |9]). In this case the generalized momen- 
tum Px — + By — is an integral of motion since 
there are no potential forces acting on electron along the 
wall (here we use the Landau gauge with a vector poten- 
tial Ax — By) . The momentum px determines a distance 
j/c between a cyclotron center and the wall, which also re- 
mains constant in time. The Hamiltonian of the system 
has the form: 

H =pI /2+ {px- By)y2 + e^y cos ujt + VUv) , (2) 

where Vw{y) is the wall potential being zero or infinity 
for y < or y > 0. Thus, we have here a so called case 
of one and half degrees of freedom (due to periodic time 
dependence of Hamiltonian on time) and the Poincare 
section has continuous invariant curves in the integrable 
regions of phase space [131 [T4L 

The Poincare sections forQ, ^ a.t j = 7/4,9/4 and 
various amplitudes of microwave field e are shown in 
Fig. [1] It shows a velocity Vy at moments of collision 
with the wall at y = as a function of microwave phase 
(f> — ujt at these moments of time. All orbits initially 
start at the wall edge y = with the initial velocity 
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FIG. 1; Poincare sections of Hamiltonian ||2| for j = 7/4 (left 
column) and j = 9/4 (right column) and different amplitudes 
of microwave field e = 0.02,0.04,0.2 (from bottom to top). 
Here, the integral px/mvp = 1, trajectories start from wall 
with fixed Vx = Px = vq = 1. Data for model (Wl) at 70 — 0, 
a« = 0. 



Vx = vq = Px — He- The value of Px = He is the inte- 
gral of motion. However, the kinetic energy of electron 
— {vl + Vy)/"^ varies with time. We see that at a 
small e = 0.02 the main part of the phase space is cov- 
ered by invariant curves corresponding to integrable dy- 
namics. However, a presence of chaotic component with 
scattered points is also visible in a vicinity of separatrix 
of resonances, especially at large e = 0.2. The points at 
Vy close to zero correspond to orbits only slightly touch- 
ing the wall, while the orbits at Vy/vo ^ 1 have a large 
cyclotron radius and collide with the wall almost per- 
pendicularly. There are also sliding orbits which have 
the center of cyclotron orbit inside the wall {yc > 0) but 
we do not discuss they here. Indeed, the orbits, which 
only slightly touch the wall {jjc ~ —vf/oJc), play the most 
important role for transport since the scattering angles 
in the bulk are small for high mobility samples and an 
exchange between bulk and edge goes via such type of 
dominant orbits [S]- 

We note that the section of Fig.[l]at j = 9/4, e = 0.02 
is in a good agreement with those shown in Fig. lb of [9]. 
However, here we have single invariant curves while in 
[S] the curves have a certain finite width. This happens 



due to the fact that in ^ the Poincare section was done 
with trajectories having different values of the integral 
Px — Dc that gave some broadening of invariant curves. 
For a fixed integral value we have no overlap between 
invariant curves as it is well seen in Fig. [l] here. 

The phase space in Fig. [l] has a characteristic reso- 
nance at a certain Vy/vQ value which position depends 
on J [5j . An approximate description of the electron dy- 
namics and phase space structure can be obtained on 
a basis of the Chirikov standard map [131 In 
this description developed in [9] an electron velocity has 
an oscillating component 5vy ~ esinoji (assuming that 
uj > ujc) and a collision with the wall gives a change 
of modulus of Vy by 2dvy (like a collision with a mov- 
ing wall). For small collision angles the time between 
collisions is At — 2{n — Vy)/ujc. Indeed, 2'it/uJc is the cy- 
clotron period. However, the time between collisions is 
slightly smaller by an amount 2vy/uJc'- at Vy <^ Vx ~ vp 
an electron moves in an effective triangular well created 
by the Lorentz force and like for a stone thrown against 
a gravitational field this gives the above reduction of At 
(formally this expression for At is valid for sliding orbits 
but for orbits slightly touching the wall we have the same 
At but with minus that gives the correction —2vy/uJc). 
The same result can be obtained via semiclassical quan- 
tization of edge states developed in jTHj. It also can be 
found from a geometric overlap between the wall and 
cyclotron circle. This yields an approximate dynamics 
description in terms of the Chirikov standard map [13j : 

Vy — Vy + 2e sin (f) + Ice, <t> — <t> + 2{-K ~ VyI p)uj/ljJc , (3) 

with the chaos parameter K — 4ea;/ (pujc)- Usually we are 
in the integrable regime with K < 1 due to small values 
of e used in experiments. A developed chaos appears 
a.t K > 1 [121 m]. Here bars mark the new values of 
variables going from one collision to a next one, Vy is the 
velocity component perpendicular to the wall, (f> = cut is 
the microwave phase at the moment of collision. Here we 
introduced a dimensionless parameter p which is equal to 
p = 1 for the case of the wall model W2 considered here. 
However, we will show that for the dynamics around disk 
with a radial field in model (DRl) we have the same map 
([3]) with p — 1 + Tc/rd- Due to that it is convenient to 
write all formula with p. We note that a similar map ^ 
describes also a particle dynamics in a one-dinicnsional 
triangular well and a monochromatic field ^17; . 

The term Ice — —IcVy + ck„ in (m describes dissipa- 
tion and noise. The later gives fluctuations of veloc- 
ity Vy at each iteration {—a < a„ < a; corresponding 
to random rotation of velocity vector in ([l])). Damp- 
ing from electron-phonon and electron-electron collisions 
contribute to 7c. The Poincare sections of this map are in 
a good agreement with those obtained from the Hamilto- 
nian dynamics as it is seen in Fig. [T] here and Fig.l in [S]. 
Following [9] we call this system model (W2) (equivalent 
to model (2) in [9]). 

A phase shift of by 27r does not affect the dynam- 
ics and due to that the phase space structure changes 
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periodically with integer values of j. Indeed, the posi- 
tion of the main resonance corresponds to a change of 
phase by an integer number of 27r values (j) — (j) — 27rm — 
2(7r — Vy/p)uj/ujc that gives the position of resonance at 
Vres — Vy — 7r/9(l — mujc/uj) — npSj/ j where m is the 
nearest integer of uj/uJc and 6j is the fractional part of 
j. Due to this relation we have the different resonance 
position for j = 7/4 and 9/4 being in agreement with the 
data of Fig. [l] at small values of e when nonlinear correc- 
tions are small (we have here p = 1). Thus at j = 9/4 
we have the resonance position at Vy = CIIIItt « 0.35 
in agreement with Fig. [I] (right bottom panel). For j = 2 
we have Vy — and at j — 7/4 the resonance position 
moves to negative value Vy = —0.45. Thus, at j = 2; 7/4 
the resonance separatrix easily moves particles out from 
the edge at Uj, < where they escape to the bulk due to 
noise. In contrast at j = 9/4 particles move along sepa- 
ratrix closer to the edge being then captured inside the 
resonance which gives synchronization of cyclotron phase 
with the microwave phase. This mechanism stabilizes the 
transmission along the edge. 




FIG. 2: (Color online) Density distribution w of electrons 
as a function of their dimensionless cyclotron center position 
yc/rc between two walls and the frequency ratio j — uj/lUc- 
The distance between specular walls is 6rc. The amplitude 
of microwave field is e = 0.1 with polarization parallel (left 
panel) and perpendicular (right panel) to the walls (see text 
for more details). Here 7o/'»J = 0.05, Oi = 0.01, Ti — l/oj. 
The variation of j — lu/uJc is obtained by changing magnetic 
field (cjc ~ B) keeping cj = const; 100 electrons are simulated 
at each j up to time = lO^/cj. Density is proportional to 
color changing from zero (black) to maximal density (white). 
Data for model (Wl). 



In [2] it is shown that the orbits started in edge vicin- 
ity are strongly affected by a microwave field that leads 
to ZRS type oscillations of transmission along the edge 
and longitudinal resistivity Rxx- The ZRS structure ap- 
pears both in the frame of dynamics described by (fl]) 
(model (Wl)) and map description ([3| (model (W2)). 
The physical mechanism is based on synchronization of a 
cyclotron phase with a phase of microwave driving that 
leads to stabilization of electron transport along the edge. 
An extensive amount of numerical data has been pre- 
sented in and we think there is no need to add more. 
Here, we simply want to illustrate that even those or- 
bits which start in the bulk are affected by this synchro- 
nization effect. For that we take a band of two walls 



with a band width between them being Aj/ = L = Qr^. 
Initially 100 trajectories are distributed randomly in a 
bulk part between walls when a cyclotron radius is not 
touching the walls (— 2rc < Uc < 2rc). Their dynam- 
ics is followed during the run time tr — lO'^/uJ accord- 
ing to Eq. ([T]) and a density distribution w{jjc) averaged 
in a time interval 5 • 10* < wt < 10^ is obtained for a 
range of 0.5 < j < 7 (261 values of j are taken homoge- 
neously in this interval). The value of tr approximately 
corresponds to a distance propagation along the wall of 
T^w ^ Vytr ^ O.lvptr ^ 5 • lO^vp/tu ^ 0.2cm at typi- 
cal values vp ^ 2 ■ lO^cm/s, Ci;/27r = lOOGHz. This is 
comparable with a usual sample size used in experiments 
[H [2] . Similar values of r^, were used in |9] . 




FIG. 3: Profile of density distribution w(yc) as a function 
of yc/rc for microwave polarization parallel (left panels) and 
perpendicular (right panels) to the walls. Here we have no 
microwave at top panels, e = 0.1, j = 1.7 at middle panels, 
e = 0.1, j — 2.4 at bottom panels. In all panels we have noise 
amplitude Qi = 0.01 as in Fig. [2] dashed curves in bottom 
panels are obtained with — 0.05. Simulations are done 
with 500 trajectories, other parameters are the same as in 
Fig. [2] Data for model (Wl). 

The dependence of density w on j/c and j are shown 
in Fig. [2] for two polarizations of microwave field. The 
data show that orbits from a bulk can be captured in 
edge vicinity for a long time giving an increase of density 
in a vicinity of edge. This capture is significant around 
resonance values j ~ jr- This is confirmed by a direct 
comparison of density profiles in Fig. [3] at j = 1.7 « 
2-1/4 and j = 2A ^ 2 + 1/4. In the later case we 
have a large density peak due to trajectories trapped in 
a resonance (see Fig. [I]) where they are synchronized with 
a microwave field. An increase of noise amplitude ai gives 
a significant reduction of the amplitude of these resonant 
peaks (Fig. [s] bottom panels). The increase of density 
is more pronounced for polarization perpendicular to the 
wall in agreement with data shown in Fig. 2 of [3]. 

We also performed numerical simulations using Eq. ^ 
with a smooth wall modeled by a potential Vw{y) — 
For large values k/uJc (e.g. k/wc = 10) we find 
the Poincare sections to be rather similar to those shown 
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m Fig. [1] that gives a similar structure of electron density 
as in Figs. |2|3[ A finite wall rigidity can produce a cer- 
tain shift of optimal capture conditions appearing as a 
result of additional correction to a cyclotron period due 
to a part of orbit inside the wall. 

The data presented in this Section show that electrons 
from the bulk part of the sample can be captured for 
a long time in edge vicinity thereby increasing the elec- 
tron density near the edge. This effect is very similar to 
the accumulation of electrons on the edges of the electron 
cloud under ZRS conditions that was reported for surface 
electrons on Helium in [7|. However we have to empha- 
size that the confinement potential for surface electrons 
is very different from the hard wall potential assumed 
in our simulations, as a consequence our results cannot 
be applied directly to this case. It is possible that the 
formation of ballistic channels on the edge of the sample 
combined with the redistribution of the electrons density 
can effectively short the bulk contribution and induce di- 
rectly a vanishing Rxx- However, it is also important 
to understand how a scattering on impurities inside the 
bulk is affected by a microwave radiation. We study this 
question in next Sections. 



III. SCATTERING ON A SINGLE DISK 
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It should be noted that resistivity properties of a reg- 
ular lattice of disk antidots in 2DEG had been studied 
experimentally dSUl] and theoretically dHlIlI]- But ef- 
fects of microwave field were not considered till present. 

In our studies we model an impurity as a rigid disk of 
fixed radius r^j = vp/uj keeping uj = const and changing 
LOc = B. In a magnetic field a cyclotron radius moves 
in a free space only due to a static rfc-electric field Edc- 
We fix the direction of Edc along x— axis and measure its 
strength by a dimensionless parameter ~ {<^vp). 
Even in absence of a microwave field a motion in a vicin- 
ity of disk in crossed static electric and magnetic fields 
of moderate strength is not so simple. The studies pre- 
sented in |22j and }23j show that dynamics in disk vicin- 
ity is described by a symplectic disk map which is rather 
similar to the map ([3]). It is characterized by a chaos 
parameter Cd — ^-nvdjlvd^c^ where Vd = vpEdc/B is the 
drift velocity; gives an amplitude of change of radial 
velocity at collision. Orbits from a vicinity of disk can 
escape for > 0.45 [22] . 

We start our analysis from the construction of the 
Poincare section in presence of microwave field at zero 
static field. To have a case with one and half degrees of 
freedom we start from a model case when a microwave 
field is directed only along radius from a disk center. The 
dynamics is described by Eq. ([T]) with a dimensionless 
microwave amplitude e. The dynamical evolution is ob- 
tained numerically by the Runge-Kutta method. At first 
we consider a case without dissipation and noise. Due to 
radial force direction the orbital momentum is an addi- 
tional integral of motion (as Px = t/c for the wall case) and 



FIG. 4: Poincare section for Hamiltonian dynamics in a disk 
vicinity in presence of radial microwave field. Left column 
panels: j = 7/4,2,9/4 at e = 0.04 (from top to bottom); 
right column panels: j = 7/4, 13/4, 9/4 at e = 0.02. Here the 
integral of orbital momentum is lo/vprd = vq/vf ~ 1, trajec- 
tories start from disk with fixed tangent velocity component 
vq — 1. Data for model (DRl). 



thus we have again 3/2 degrees of freedom. We call this 
disk model with radial microwave field as model (DRl). 

The Poincare sections at the moments of collisions with 
disk are shown in Fig. |4j [5] for model (DRl). Here, 
is the radial component of electron velocity and is a 
microwave phase both taken at the moment of collision 
with disk. We see that the phase space structure re- 
mains approximately the same when j is increased by 
unity (compare j ~ 9/4,13/4 panels in Fig. |4]). This 
happens for orbits only slightly touching the disk (small 
Vr) since the microwave phase change during a cyclotron 
period is shifted by an integer amount of 27r (in a first 
approximation at <C Tc). The similarity between the 
wall and disk cases is directly seen from Fig. [5] as well as 
periodicity with j — > j + I. 

In fact in the case of disk with a radial field the dynam- 
ics can be also described by the Chirikov standard map 
([3]) where Vy should be understood as a radial velocity 
Vr at the moment of collision. The second equation has 
the same form since the change of the phase between two 
collisions is given by the same equation but with the pa- 
rameter p = 1 + Tc/rd- This expression for p is obtained 
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from the geometry of slightly intersecting circles of ra- 
dius Td for disk and radius Tc for cyclotron orbit (the 
angle segment of cyclotron circle is A(/3 = 2vr/p)- For 
fd ^ "fc this expression naturally reproduces the wall 
case while at ^ Tc we have the correction term pro- 
portional to Vy going to zero that also well corresponds 
to the geometry of two disks. After such modification of 
p we find that the resonance positions Vres = T^pSj /j are 
proportional to p. Thus the model (DRl) reduced to the 
map description ([3| at p > 1 is called model (DR2). 

The expression for Vres works rather well. Indeed, for 
j = 2.1 in Fig. [5] we obtain w^es = 0.149 for model (W2) 
and 0.463 for model (DR2). These values are in a good 
agreement with numerical values Vres ~ 0.15 for model 
(W2) and Vres ~ 0.6 for model (DRl). In the later case 
the agreement is less accurate due to a larger size of non- 
linear resonance. The comparison of Poincare sections 
given by the Chirikov standard map ([s]) and the dynam- 
ics from Newton equations, shown in Fig.[5j confirms the 
validity of map description. 

According to the well established results for the 
Chirikov standard map [T3j we find for models (Wl), 
(W2) and (DRl), (DR2) the width of separatrix 5v and 
the corresponding resonance energy width E.^ = ((5v)^/2: 

Er = IGeujcpEp/u! ; p = 1 + rc/r^ ; (4) 

Vres = TTp6j/j ; 6v ^ '^\J (P 1 3 
Sje = Svj/{2TTp) ; j = 1+Uj/uJc , 

where Sje is the resonance shift produced by a resonance 
half width 6v/2 = Vres- This relation shows that for the 
disk case this energy is increased by a factor p compared 
to the wall case. In majority of our numerical simulations 
we have p = 1 + j- 

Thus a radial field models (DRl), (DR2) represent a 
useful approximation to understand the properties of dy- 
namics in a disk vicinity but a real situation corresponds 
to a linear microwave polarization and the Poincare sec- 
tion analysis should be modified to understand the dy- 
namics in this case. 

Due to that we start to analyze the scattering problem 
on a disk in presence of weak static field eg and microwave 
field e using Eq. Q. For the scattering problem we find 
more simple to have dissipation to work only at the time 
moments of electron collisions with disk: at such time 
moments the radial component of electron velocity is re- 
duced by a factor Vr — )■ Vr/(1 + 7(i), the reduction is 
done only if the kinetic energy of electron is larger than 
the Fermi energy. Such a dissipation can be induced by 
phonon excitations inside the antidot disk. We fix geom- 
etry directing rfc-field along a;— axis and microwave along 
?/— axis. The noise is modeled in the same way as above 
in Eq. ([T]). We call this system disk model (Dl). 

Examples of electron cyclotron trajectories scattering 
on disk are shown in Fig. [6] In absence of microwave 
field a trajectory escapes from disk rather rapidly. A 
similar situation appears at j = 2 and microwave field 
with e = 0.04. In contrast for j = 9/4 and e = 0.04 a 




FIG. 5: Poincare sections for wall model (Wl) (left) and disk 
model with radial electric field (DRl) (right) at e = 0.01 and 
j — 2.1 (top, middle) and j = 3.1 (bottom). Top panels 
are obtained from the Chirikov standard map (|3| at p = 1 
(left top panel), corresponding to the wall model (W2), and 
at p = 1 -|- ji = 3.1 (right top panel), corresponding to the 
disk model (DR2). Middle and bottom panels are obtained 
from solution of Newton equations |T| for wall (left) and disk 
(right). Other parameters are as in Fig. [l] and Fig. |4] Vy and 
Vr are expressed in units of vp. There is no dissipation and 
no noise. 

trajectory can be captured for a long time or even forever 
depending on initial impact parameter. 

For some impact parameters a trajectory can be cap- 
tured for a very long time tc, in certain cases in absence 
of noise we have tc — oo. At such long capture times the 
collisions with disk become synchronized with the phase 
(j) of microwave field at the moment of collisions. This 
is directly illustrated in Fig. [7] where we show the angle 
of a collision point on disk, counted from axis, in 
dependence on (j). Indeed, the dependence 6 on (j) forms 
a smooth curve corresponding to synchronization of two 
phases. At the same time the radial velocity at collisions 
Vr moves along some smooth invariant curve Vr{<j)) in the 
phase space (w^, (/>). However, to make a correct compari- 
son with the radial field models (DRl), (DR2) we should 
take into account that the cyclotron circle rotates around 
disk so that we should draw the Poincare section in the 
rotational phase (j)' ~ cj) — 6. In this representation we 
see the appearance of the resonance (see right column of 
Fig. [t]) that is similar to those seen in Figs. |4] [5] for the 
radial field models. 

A more direct correspondence between radial field 
models (DRl), (DR2) and the model (Dl) with a linearly 
polarized microwave field is well seen from the Poincare 
sections shown in the rotation frame of phase (j)' = (f) — 9 
in Fig. [8j In this frame we see directly the resonance 
at j — 2.1,2.25 being very similar to the wall case and 
the radial field model. However, the positions of reso- 
nance at Vr = Vres are different from those in Fig. [5] 
Of course in the rotation frame the orbital momentum is 
only approximately conserved that gives a broadening of 
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FIG. 6: (Color online) Scattering of electron cyclotron tra- 
jectory on a disk scatterer (blue/black circle) in model (Dl). 
Top left panel shows the case in absence of microwave, e = 0, 
j = 9/4. Top right panel: temporary captured path at 
e = 0.04, j — 9/4. Bottom left panel: path captured for- 
ever at e = 0.04, j — 9/4. Bottom right panel: no capture 
at e = 0.04, j = 2. The trajectory part colliding with disk 
is shown by black curve, its part before and after collisions is 
shown in gray. The red (light gray) points and curves show 
the trajectory of cyclotron center. Here the dissipation pa- 
rameter is 7d/cj — 0.01; the static electric field is directed 
along x— axis and = 0.001; microwave field is directed along 
t/— axis. There is no noise here. Coordinates x, y are expressed 
in units of rd- Data for model (Dl). 



invariant curves in Fig. |8] 

We explain this as follows. For the linear polarized 
field of model (Dl) the radial component of microwave 
field is proportional to ^ e sin 9 cos ^ 0.5e sin(aji— 0) 
where we kept only slow frequency component of radial 
field (the neglected term with sin((u;i -I- 9) gives resonant 
values Vres > vf)- The radial field gives kicks to the 
radial velocity component at collisions with disk similar 
to the case of model (DR2) described by Eq. ([3|: = 
Vr + 0.5esin(/)', 0' = (/>' (27rj - 2vrj/p) - 2vr{p - l)/p. 
Here we use the radial field component phase (j)' = ujt~9 
at a moment of collision with disk (the tangent compo- 
nent does not change and can be neglected). The 
phase variation 0' — (j)' has the first term 2nj — 2vrj / P 
being the same as for the radial field model (DR2), and 
an additional term related to rotation around disk with 
— = —2vr{p — l)/p which comes from geometry. In- 
deed, the segment angles of intersections of circles and 
Tc are: for disk radius rd it is A^? = 2vr{p — 1) / p and for 
cyclotron radius r,. it is A(p = 2vr/ p. Thus, their ratio is 
A0/A(^ = r^/r^ in agreement with the geometrical scal- 
ing. This result can be obtained from the expression for 




FIC. 7: Synchronization between disk collision angle 6 and 
microwave phase cf). Left column: Dependence of angle 6 of 
collision point on disk, counted from i— axis, and radial ve- 
locity Vr, taken at collision, on microwave phase (f). Right col- 
umn: Same as in left column with (p' = (j) — 9. Here j — 2.25, 
e = 0.04, = 0.001, jd/(^ ~ 0.01, vo — vp, there is no noise; 
points are shown for times W^/uj < t < lO^/o;, the capture 
time of this orbit is tc > lO^/oj. Data for model (Dl). 



0.8- 



0.4- 



0.0 

0.8- 



1 * 



"T" 




0.4- 



0.0 
0.8- 



04 



0.0- 



-1 



1 -1 




FIC. 8: Poincare sections in phase plane {vr,cf>') with cf>' = 
4>~e ior j = 2.1, e = 0.01 (left top); j = 2.25, e = 0.01 (right 
top); j = 2.75, e = 0.02 (left middle); j = 2.25, e = 0.02 
(right middle); j = 2.75, e = 0.04 (left bottom); j = 2.25, 
e = 0.04 (right bottom); Vr is expressed in units of vp- Data 
for model (Dl), no noise no dissipation. 



A(p by interchange of two disks that gives the above ex- 
pression for A9 (at = both shifts Aip = 2vr/p and 
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= 2vr{p— l)/p are equal). 
Thus again the dynamical description is reduced down 
to the Chirikov standard map with slightly modified pa- 
rameters giving us for the model (Dl) the chaos parame- 
ter K — 2e(j + p—l)/p being usually smaller than unity, 
resonance position Vres, resonance width 6v and the res- 
onance energy width Er = {5vY /2: 



SepEp/ip + J - 1); P=l + rc/rd ; (5) 
T^pSj/U + P-1); 5v^ 4^ep/{2{j + p-l)); 
6v{p + j - l)/(27rp); j = uj/uJc , 



where Sj^ is a shift of resonance produced by a finite 
separatrix half width 6v/2. For our numerical simula- 
tions we have p = 1 + j with Vres = T^{j + l)'5j/(2j), 
5v = 2^e{] + and 5j, = (2/^) ^^jVO + 1). 

At p = J + 1 Eq. ([5]) gives the values w^es = 0.232 at 
j = 2.1 while the numerical data of Fig.[8]give Vres ~ 0.2, 
and we have at j = 2.25 the theory value Vres = 0.567 
being in a good agreement with the numerical value 
Vres ~ 0.5 of Fig. [8j For j — 2.75 we have the reso- 
nance position at Vr < corresponding to the bulk and 
thus the resonance is absent. The resonance width in 
Fig. [8] at j = 2.25, e = 0.01 can be estimated as Sv ~ 0.3 
that is in a satisfactory agreement with the theoretical 
value Sv = 0.24 from ([s]). We remind that in model 
(Dl) we have only approximate conservation of orbital 
momentum that gives a broadening of invariant curves 
and makes determination of the resonance width less ac- 
curate. In spite of this broadening we see that the res- 
onance description by the Chirikov standard map works 
rather well. 

In Fig. [9] we show the dependence of average capture 
time tc on j in the model (Dl). The averaging is done 
over Ng = 500 trajectories scattered on disk at all such 
impact parameters that cyclotron orbit can touch the 
disk. Here, we show the ratio tc/tc{0) where tc{0) is an 
average capture time in absence of microwave. According 
to our numerical data we have approximate dependence 
ujtc{0) ~ 3/y^ corresponding to a period of nonlinear 
oscillations in a disk map description discussed in j22[|23j . 

The data of Fig. [9] show a clear periodic dependence 
of capture time tc on j corresponding to the periodicity 
variation of Poincare section with j (see Figs. [4j [5] [s]). 
This structure is especially well visible at weak static 
fields. With an increase of €s this structure is suppressed. 
Indeed, at large even without microwave field the tra- 
jectories can escape from disk as it discussed in j^H US] 
and microwave field does not affect the scattering in this 
regime. 

The distributions of capture times are shown in Fig.fTol 
We clearly see that at resonant values of j a microwave 
field leads to appearance of long capture times. For 
example, we have the probability to be captured for 
tc > 180/w being W = 0.46 at j = 2.25 while at j = 2 we 
have W = 0(left panel in Fig.[w]); and we have W = 0.38 
at j = 2.37 while at j = 1.9 we have W <3- 10"^ (right 




FIG. 9: Dependence of rescaled capture time tc/tc{0) on j 
at e = 0.04 shown at various static fields: £b = 2.5 • 10"'* (top 
left), 10"^ (top right), 8- 10"^ (bottom left), 0.016 (bottom 
right). Here 7^/1^; = 0.01, there is no noise. Data for model 
(Dl). 



panel in Fig. 10 1 . These data confirm much stronger cap- 



ture at certain resonant values of j. 




FIG. 10: Differential probability distribution t^dW/dtc of 
capture times tc for parameters of Fig. |9] Left panel: j = 
2 (dashed curve for minimum of capturing probability) and 
j — 2.25 (full curve at maximum of capturing probability) 
at 65 — 0.001. Right panel: similar cases at j = 1.9 (dashed 
curve) and j — 2.37 (full curve) at = 0.002. Data are 
obtained with Ns = 5 ■ 10'* trajectories started at different 
impact parameters and running up to time t = 3 ■ 10^ /cj. 
Here tc is expressed in units of 1/uj. Data for model (Dl). 



The data of Figs. |9j [TO] show that the scattering pro- 
cess on disk is strongly modified by a microwave field. 
However, to determine the conductivity properties of a 
sample we need to know what is an average displacement 
Ax along static field after a scattering on a single disk. 
Indeed, in our model a dissipation is present only during 
collisions with disk while in a free space between disks the 
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dynamics is integrable and Hamiltonian. Hence during 
such a free space motion there is no displacement along 
the static field (the dissipative part of conductivity or 
resistivity appears only due to dissipation on disk). The 
dependence of Ax on j is shown in Fig. In absence 
of microwave field at e = we find Ax oc 1/j cx _B that 
corresponds to a simple estimate Ax cx LOdd- The nu- 
merical data show that Ax is practically independent of 
Es and that Ax = in absence of dissipation at 7^ = 0. 
In presence of microwave field we see that the displace- 
ment along the static field has strong periodic oscillations 
with j. The striking feature of Fig. [TT]is the appearance 
of windows of zero displacement Ax ~ at resonance 
values jr = 9/4, 13/4, 17/4.... We discuss how this scat- 
tering on a single disk modifies the resistance of a sample 
with large number of disks in next Section. 
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FIG. 11: Average shift Ax along static field after scattering 
on a single disk shown as a function of j at various amplitudes 
of static and microwave fields: = 0.0005, e = (gray dashed 
curve); ts = 0.001, e = (black dashed curve); = 0.0005, 
e = 0.04 (gray fuU curve); = 0.001, e = 0.04 (black full 
curve). The data are obtained by averaging over Ns — 5 ■ 10^ 
scattered trajectories with random impact parameters; here 
7d/a; — 0.01, noise amplitude — 0.005. Data for model 
(Dl). 



IV. RESISTANCE OF SAMPLES WITH MANY 
DISKS 

To determine a resistance of a sample with many disks 
we use the following scattering disk model. The scatter- 
ing on a single disk in a static electric field £5 is com- 
puted as it is described in the previous Section with a 
random impact parameter inside the collision cross sec- 
tion ad = 2{rc+rii). After that a trajectory evolves along 
y-axis according to the exact solution of Hamiltonian 
Eq. ([I]) (no dissipation and no noise) up to a collision with 
next disk which is taken randomly on a distance between 
2(rc -I- Trf) and 21^ where = i/{o'dnd) is a mean free 
path along y— axis and Ud is a two-dimensional density of 



disks (of course £e ^ 2{rc + rd)). In a vicinity of disk the 
dynamical evolution is obtained by Runge-Kutta solution 
of dynamical equations as it was the case in previous Sec- 
tion. We use low disk density with Udr'^ ^ 1/100. The 
collision with disk is done with a random impact param- 
eter in the x— axis of disk vicinity: the impact parameter 
is taken randomly in the interval [— (r^ -I- rd), {r^ + rd)] 
around disk center. Noise acts only when a center of cy- 
clotron radius of trajectory is on a distance r < rd + Tc 
from disk center so that a collision with disk is possible. 
After scattering on a disk a free propagation follows up 
to next collision with disk. 

Along such a trajectory we compute the average dis- 
placement 5x and 6y after a time interval 5t. In this 
way the number of collisions with disks is Ncoi ~ Sy/£e 
and a total displacement in x— axis is 6x ~ N^-Ax where 
Ax is an average displacement on one disk discussed in 
previous Section. We compute the global displacements 
6x, Sy on a time interval 6t — 10® /a; averaging data over 
200 trajectories. We call this system model (D2). 
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FIG. 12: Dependence of resistivity Rxx and R^y on magnetic 
field B — 1/j (resistivity is expressed in arbitrary numeri- 
cal units). Blue and gray curves show respectively Rxx and 
Rxy in absence of microwave field. Black curve with points 
show Rxx dependence of B at microwave field e = 0.04 (here 
B is expressed in units of l/o; = const). Here es = 0.001, 
7ci/ai = 0.01, noise amplitude Oi = 0.005, = l/uj. Res- 
onant values jr are shown by arrows; bars show statistical 
errors for Rxx- Data are obtained by averaging over 200 tra- 
jectories propagating up to time t = 10^ /cu. Data for model 
(D2). 

Then the current components are equal to jx — 5x/St, 
jy = Sy/St and conductivity components are axx — 
jx/Edc, a-xy = jy/Edc (the current is computed per 
one electron). We work in the regime of weak rfc-field 
where jxijy scales linearly with Edc- The current jy 
is determined by the drift velocity Vd = Edc/B ^ vp- 
Since the mean free path is large compared to disk 
size 3> > rd we have an approximate relation 



'xy 



1/B, 



Ax/B£e « axyA. 



As 



2DEG experiments [H |2] we have in our simulations 



'xy 



/o-x 



100 (see Fig. 12). The resis- 



tivity is obtained by the usual inversion of conductiv- 
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ity tensor with 



Rxy = <^xyl{<yxx + ^xy) ~ '^/<^xy The dependence oiRxx, 
Rxy, expressed in arbitrary numerical units, on magnetic 
field B = u}/j = 1/j is shown in Fig. [12} 

In absence of microwave field we find Rxy oc B and 
Rxy/Rxx ~ 200 similar to experiments For small 

noise amplitude (e.g. a; — 0.005) we have Rxx growing 



linearly with B (see Fig. 12 ) but at larger amplitudes (e.g. 
at ~ 0.02) its increase with B becomes practically flat 
showing only 30% increase in a give range of B variation. 
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FIG. 14: (Color online) Same as in Fig. [13] at various 
microwave amplitudes e — 0.01 (red/gray curve), 0.02 
(blue/dark curve), 0.04 (black curve); amplitude of noise is 
fixed at at = 0.005. Data for model (D2). 



tion of cyclotron phase with microwave field phase and 
capture inside the nonlinear resonance. 



FIG. 13: (Color online) Rescaled values of resistivity R^x 
(left panel) and R^y (right panel) as function of j = uj/oJc 
at various noise amplitudes — 0.005 (black curve), 0.01 
(blue/dark curve), 0.02 (red/gray curve). Here e = 0.04, = 
0.001, other parameters are as in Fig. |12[ Curves are drown 
though numerical points obtained with a step Aj = 0.1. Data 
for model (D2). 

In presence of microwave field the dependence of R^x 
on B is characterized by periodic oscillations with mini- 
mal Rxx values being close to zero at resonant values of 
j = jr well visible in Fig. [12] The dependence of Rxx , 
Rxy rescaled to their values Rx 



(0), Rxy{0) in absence 



of microwave field are shown in Fig. |13|a t various ampli- 
tudes of noise and fixed e, and in Fig. |14| at various e and 
fixed noise amplitude a^. We see that increase on noise 
leads to an increase of minimal values of Rxx at resonant 
values jr- In a similar way a decrease of microwave power 
leads to increase of minimal values of Rxx at jr- At the 
same time the Hall resistance Rxy is only weakly affected 
by microwave radiation as it also happens in ZRS exper- 
iments. 

These results are in a qualitative agreement with the 
ZRS experiments. On the basis of our numerical stud- 
ies we attribute the appearance of approximately zero 
resistance at jr values in our bulk model of disk scat- 
terers to long capture times of orbits in disk vicinity at 
these jr values (see Fig. [9]). During this time tc noise 
gives fluctuations of coUisional phase 9 and due to that 
a cyclotron circle escapes from disk practically at ran- 
dom displacement Aa; that after averaging gives average 
Aa; = 0. Since resistivity is determined by the average 
value of Ax this leads to appearance of ZRS. We note 
that this mechanism is different from the one of edge 
transport stabilization discussed here and in [9]. How- 
ever, both mechanisms are related to a long capture times 
near edge or near disk that happens due to synchroniza- 
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FIG. 15: (Color online) Phase space {vr, (j)') of trajectories at 
the moment of collisions with disks for parameters of Fig. |l4] 
at e = 0.04: j = 2.1 (top left), j = 2.25 (top right), j = 2.75 
(bottom left), j = 3.25 (bottom right). Here 0' = <j) — 9 where 
(j) = u)t is a, microwave phase at the moments of collisions 
with disk and 9 is the angle on disk at collision moment, 
counted from x-axis (same as in Figs.[7][8|. Data are obtained 
from 500 trajectories iterated up to time t = lO^/u;. Density 
of points is shown by color with black at zero and white at 
maximum density. The average number of collisions per disk 
per trajectory is N^oi = 12.4, 25.9, 9.5, 15.5 respectively for 
j = 2.1, 2.25, 2.75, 3.75. Data are obtained for model (D2). 

To illustrate the capture inside the resonance we 
present the distributions of trajectories from Fig. [T4| 
shown in the phase space plane (vr^cj)') at the moments 
of collisions with disks in Fi g. \TE\ This is similar to 
the Poincare sections of Fig. [8| however, now we con- 
sider the real case of diffusion and scattering on many 
disks in the model (D2) with noise and dissipation. We 
see that for j = 2.25 orbits are captured in a vicinity 
of the center of nonlinear resonance at 0' well seen 
m Fig. [8} For j = 2.1 we have a density maximum lo- 
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cated at smaller values of Vres and 



even if there 



is a certain shift of Vres produced by a significant reso- 
nance width at e = 0.04. At = 2.75 we have a density 
maximum at 0' « ±7r corresponding to an unstable fixed 
point of separatrix. The total number of collision points 
Ncoi in this case is by a factor 2.5 smaller than in the 
case of stable fixed point at j — 2.25. A similar situation 
is seen in the case of wall model (Wl) (see Fig. Id, f in [9 ) 
even if there the ratio between number of captured points 
was significantly larger. The results of Fig. [15] show that 
in the ZRS phase the collisions with disk indeed create 
synchronization of cyclotron and microwave phases and 
capture of trajectories inside the nonlinear resonance. 

However, there are also some distinctions between bulk 
disk model (D2) and experimental observations. The first 
one is that there are minima for Rxx/ Rxxi^) but there 
are no peaks which are well visible around integer j val- 
ues in ZRS experiments [B [2] i [S] and numerical simula- 
tions of transport along the edge [9]. The second one is 
appearance of small negative values of Rxx at jV values. 
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FIG. 16: Dependence of Rxx/RxxiQ) on j in the disk model 
with dissipation at disk collisions at rate 'yd/^ = 0.01 and dis- 
sipation in disk vicinity with rate 'yo/uj = 0.02; a disk rough- 
ness gives additional angle rotations with amplitude ad = 0.1 
(see text); the amplitude of noise in disk vicinity is ai = 0.001. 
Here we have e = 0.04, Es — 0.001; 51 numerical points in j 
are connected by lines to adapt an eye. Data are obtained 
by averaging over 100 trajectories propagating up to time 
t = W^/lo. Data for model (D3). 

We attribute the absence of peaks to a specific dissipa- 
tion mechanism which takes place only at disk collisions. 
It is rather convenient to run long trajectories using exact 
solution for free propagation between disks. Indeed, in 
this scheme there is no dissipation during this free space 
propagation and thus these parts of trajectories have no 
displacement along static field. We also tested a dissi- 
pation model with additional j{v) = Jo{\v/vf\'^ — 1) for 
\v\ > vp and ^(v) = for \v\ < vp. This dissipation 
works only in a disk vicinity when the distance between 
disk center and cyclotron center is smaller than + rc. 



The dissipation 7^ on disk remains unchanged. We also 
added a certain roughness of disk surface modeled as an 
additional random angle rotation of velocity vector in the 
range ±ad, done at the moment of collision with disk. 
We call this system disk model (D3). The results for the 
resistivity ratio Rxx/ Rxx{0) are shown in Fig. |16[ They 
show an appearance of clear peaks of Rxx/Rxxi^) in pres- 
ence of such additional dissipation in vicinity of integer 
j. There is also a small shift of minima from integer plus 
1/4 to integer plus 0.4 — 0.5. 

The second point of distinctions from ZRS experiments 
is a small negative value of Rxx at resonant j values. It 
is relatively small for disk model (D2) (see Figs. [13 14) 
and it becomes more pronounced for disk model (D3) 



(see Fig. 16). It is possible that a scattering on disk in 
presence of dissipation, noise, static and microwave field 
gives a negative displacement Ax which generates such 
negative Rxx values. We expect that in the limit of static 
field going to zero this effect disappears. Indeed, the 
negative values become smaller at smaller £5 according 
to data of Fig. 11 but unfortunately the small limit is 
also very difficult to investigate numerically. 

We consider that at this stage of the theory the pres- 
ence of negatives values for Rxx does not constitute a 
critical disagreement. The escape parameters for elec- 
trons that have been captured on an impurity for a time 
long enough to make many rotations around it, are likely 
to strongly depend on the model for the electron impu- 
rity interaction and further theoretical work on a more 
microscopic model is needed. In general a zero average 
displacement along the field direction seems natural for 
a smooth distribution of trapping times with a charac- 
teristic time scale much larger than the rotation time 
around the impurity (this assumption does not seem to 
hold for our model, see for example the sharp features on 
Fig. 10). Finally in Section VII we propose a slightly dif- 
ferent mechanism by which the combination of trapping 
on impurities investigated here and electron-electron in- 
teractions can lead to ZRS. 



V. PHYSICAL SCALES OF ZRS EFFECT 

The ZRS experiments [TJ [5] show that the resistance 
Rxx in the ZRS minima scales according to Arrehenius 
law Rxx (X exp(— Tq/T) with a certain energy scale de- 
pendent on a strength of microwave field. In typical ex- 
perimental conditions one finds very large Tg w 20K at 
jr = 5/4 (see e.g. Fig. 3 in l^). These data also indicate 
the dependence Tq cx 1/jV oc B. This energy scale ksTg 
is very large being only by a factor 7 smaller than the 
Fermi energy Ep/ks ~ 150i^r. At the same time the am- 
plitude of microwave field is rather weak corresponding to 
e « 0.003 at field of IV/ cm or ten times larger at lOV/cm 
(unfortunately it is not known what is an amplitude of 
microwave field acting on an electron). 

As in [S] we argue that the Arrehenius scale is de- 
termined by the energy resonance width Q with Tq = 
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Er/ks- Indeed, the resonance forms an energy barrier 
for a particle trapped inside the resonance by dissipative 
effects being analogous to a wash-board potential. An 
escape from this potential well requires to overcome the 
energy leading to the Arrehenius law for R^x depen- 
dence on temperature. Assuming the case of the wall 
with p = 1 we obtain at E = 3V/cm the activation 
temperature Tq w 23K being is a satisfactory agreement 
with the experimental observations. The theoretical re- 
lation Q also reproduces the experimental dependence 
To (X ifjr at p = 1. In his relation Tq oc e (x E be- 
ing confirmed by the numerical simulations presented in 
[H]. This dependence is in a satisfactory agreement with 
the power dependence found in experiments [1] . In other 
samples one finds that the dependence Tq cx works in 
a better way. We think that higher terms in a nonlinear 
resonance can be responsible for scaling Tq cx being 
different from the relation Q. Also a finite rigidity of 
the wall or disk scatterers can be responsible for appear- 
ance of higher power of e. 

The energy scale Er on disks is enhanced by a factor 
p — 1 + Tc/rd for the case of radial field Q. However, 
we showed that for a linear polarization the scale Er is 
given by Eq.([5| and thus there is no enhancement at 
large p. Indeed, we performed direct simulations at pa- 
rameters of Fig. [14] with the reduced value of disk radius 
by a factor 2. The numerical data give approximately 
the same traces Rxx/Rxx{0) vs. j at e = 0.01,0.02,0.04 
without visible signs of deeper minima at small e. This 
confirms the theoretical expressions ([5). In any case, for 
small values <ti Vc one should analyze the quantum 
scattering problem which is significantly more complex 
compared to the classical case. We may assume that 
in a quantum case one should replace by a magnetic 
length aB ~ \fv oc \fB. In such a case we are getting 
p = 1 + ^i; « 9 that gives Tq « %K at j « 2.25 and 
microwave amplitude E « "iV jcm. However, in this case 
we obtain the scale Tq being practically independent of 
j which differs from experimental data. In any case in 
experiments the size of impurities is small compared to 
Tc and a quantum treatment is required to reproduce the 
correct picture for Kxx dependence on parameters in the 
ZRS phase. 

Another point is related to the positions of ZRS min- 
ima on j axis. We remind that that for the wall model the 
resonance is located at Vres = T^^j/j Q and that the sep- 
aratrix width is 5v — '^i^/Tfj. The capture of trajectories 
from the bulk is most efhcient when a half width of sepa- 
ratrix touch the border of bulk at = with v^es = ^'^ /2 
that gives the expression 5je = l^ft] /tt for the wall case. 
At e = 0.06, i = 2.25 this gives hj^ = 0.22 being in a 
good agreement with the numerical data 8j 1/4 for 
Rxx dependence on j (see Figs. 2, 3 in [9 with a visible 
tendency of growth with j). For the data presented 
here in Fig. [14] for the disk case at e = 0.04, p = 1 -I- j, 
j = 2.25 we obtain from ([5| (5je = 0.11 that is slightly less 
than the numerical value Oj' w 0.25 for minima location. 
We attribute this difference to an approximate nature of 



expression for the resonance width at relatively strong 
microwave fields. We also note that in experiments an 
additional contribution to the value of &j can appear due 
to a finite rigidity of disk and wall potentials. 



VI. THEORETICAL PREDICTIONS FOR ZRS 
EXPERIMENTS 

The theoretical models presented here and in |9j re- 
produce main experimental features of ZRS experiments 
[TJ[5],[5j. However, it would be useful to have some ad- 
ditional theoretical predictions which can be tested ex- 
perimentally. A certain characteristic feature of both 
wall and disk models is appearance of nonlinear reso- 
nance. For example, according to the wall model (W2) 
described by Eq. ([s]) the dynamics inside the resonance is 
very similar to dynamics of a pendulum. The frequency 
of phase oscillations inside the resonance is fip^ 



K : 



2yjeu]/u]c = 2^/1] <C 1 irS]. Here the frequency is ex- 
pressed in number of map iterations and since the time 
between collisions is approximately jixJc we obtain the 
physical frequency fi^ of these resonant oscillations be- 
ing r^r/w = f2p;ia;c/(27rti;) = ^ eWc/w/tt. At e ~ 0.02 this 
frequency is significantly smaller than the driving mi- 
crowave frequency. The dynamics inside the resonance 
should be very sensitive to perturbations at frequency 
wi K, that gives: 



(6) 



To check this theoretical expectation we study numer- 
ically the effect of additional microwave driving with di- 
mensional amplitude ei (ei ^ e) and frequency ui. We 
use the wall model (W2) based on the Chirikov standard 
map described here and in [9 . As for the additional driv- 
ing frequency we have ei — \Ei\/ {uivp) where Ei is the 
field strength of microwave frequency uji , we assume that 
both main and additional fields are collinear and perpen- 
dicular to the wall. In presence of second frequency the 
map ([3| takes the form: 

Vy — Uj, + 2e sin -|- 2ei sin(wi0/a;) -|- /cc) (7) 

4> — <f) + 2(tI — Vy)ljJ /UJc ■ 

The only modification appears in the first equation since 
now the change of velocity at collision depends on both 
fields; the second equation remains the same as in (|3|. 
As in the model (W2) the term Ice describes the effects 
of dissipation with rate 7c and noise with amplitude a 
of random velocity angle rotations. We call this system 
model (W3). 

In the model (W3) the resistance Rxx is computed 
numerically in the same way as in the model (W2) de- 
scribed in [9 : the displacement along the edge between 
collisions is 5x = 2vy/ojc] it determines the total dis- 
placement Ax along the edge during the total computa- 
tion time At ~ 10" /w; then Rxx oc l/Dx = At/(Aa;)^ 
where Dx is an effective diffusion rate along the edge. 
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FIG. 17: (Color online) Left panel: Dependence of rescaled re- 
sistance Rxx / RxxiO) on frequency ratio uji/uu in model (W3) 
described by the map ([7|. Here, the test driving at frequency 
uji has fixed amplitude ei = 0.007; the main microwave driv- 
ing is located in the ZRS phase at j — uj/uJc = 2.25, and 
its amplitude takes values e = 0.02 (blue curve), 0.03 (green 
curve), 0.04 (magenta curve), 0.06 (red curve) (these curves 
follow from top to bottom at cji/oj — 0.2). The values of 
Rxx are computed at fixed ei — 0.007 and corresponding e; 
the values of Rl.^{0) are computed at ei = and e = 0.02. 
The data are obtained at noise amplitude a — 0.02 and dis- 
sipation 7c = 0.01; averaging is done over 2000 trajectories 
for 5000 iterations of map ([7|. Right panel: Dependence of 
peak position uji/lj on main microwave driving amplitude e 
obtained from data of left panel at ei — 0.007 and additional 
data at ei — 0.003 (blue points), the theory dependence Q 
at lo/loc ~ 2.25 is shown by the straight red line. Data for 
model (W3). 



edges. If the dominant contribution is given by bulk disk 
scatterers then a certain numerical coefficient A should 
be introduced in the right part of the expression. Ac- 
cording to the data of Fig. |8] and Eqs. ([s]) we estimate 
A « 0.5 (the separatrix width is smaller for the disk case 
compared to the wall case at the same e). 

Another interesting experimental possibility of our the- 
ory verification is to take a Hall bar of high mobifity 
2DEG sample and put on it antidots with regular or dis- 
ordered distribution (it is important to have no direct 
collision-less path for a cyclotron radius in crossed dc- 
electric and magnetic fields) with a low density of anti- 
dot disks ridr^ <C 1 (as in our numerical studies) so that 
an average distance between antidots is larger than the 
cyclotron radius re- The regular antidot lattices have 
been already realized experimentally [TSj [19] . The effect 
of microwave field on electron transport in a regular lat- 
tice has been studied in the frame of ratchet transport 
in asymmetric lattices [25) . Even a case of symmetric 
circular antidots has been studied in [25] but the lattice 
was regular and no special attention was paid to analysis 
of resistivity at ZRS resonant regime with j ^ jr- We 
think that the experimental condition of |25j can be rel- 
atively easy modified to observe the ZRS effect on disk 
scatterers discussed here. 



VII. DISCUSSION 



To see the effect of additional weak test driving ei at 
frequency uji we place the system in the ZRS phase at 
j = uj/uJc = 2.25 and measure the variation of rescaled re- 
sistance Rxx/Rxxi^)- Here Rxx is the resistance in pres- 
ence of both microwave fields e and ei while Rxxi^) is 
the resistance at ei = and a certain fixed e. The de- 
pendence of Rxx/R%x{^) on the frequency ratio is shown 



in Fig. 17 at left panel. The main feature of this data is 
appearance of peak at low frequency ratio wi/w < 0.1. 
In the range 0.1 < wi/oj < 0.4 the testing field ei is non- 
resonant and does not affect Rxx however at uji/lu < 0.1 
it becomes resonant to the pendulum oscillations in the 
wall vicinity and hence strongly modifies Rxx value. The 
dependence of this resonance ratio wi /ui on amplitude of 



main driving field e is shown in right panel of Fig. 17 
The numerical data are in a good agreement with the 
above theoretical expression m. 

The theoretical dependence^6]) allows to check the syn- 
chronization theory of edge state stabilization. It also al- 
lows to measure the strength of main microwave driving 
force acting on an electron that still remains an exper- 
imental challenge. The experimental testing of relation 
([g]) requires to work with good ZRS samples which have 
every low resistance in ZRS minima since this makes the 
effect of testing field wi to be more visible. We note 
that the recent experiments in a low frequency regime 
uj/uJc 1 [24] demonstrate that Rxx is sensitive to low 
frequency driving. The expression ([6]) is written for the 
case when Rxx is mainly determined by a transport along 



Above we presented theoretical and numerical results 
which in our opinion explain the appearance of microwave 
induced ZRS in high mobility samples. The synchroniza- 
tion theory of ZRS proposed in [3] and extended here 
is based on a clear physical picture: high harmonics 
w/Wc = J > 1 are generated by collisions with sharp edge 
boundary or isolated impurities which are modeled here 
by specular disks. The ZRS phases appear in a vicinity 
of resonant values j,- ~ 1 + 1/4, 2-1- 1/4, .... At these 
ir values the cyclotron phase of electron motion becomes 
synchronized with the microwave phase due to dissipative 
processes present in the system. 

For trajectories at the edge vicinity this synchroniza- 
tion gives stabilization of propagation along edge chan- 
nels that creates an exponential drop of resistivity con- 
tribution of these channels with decreasing amplitude 
of thermal noise and increasing amplitude of microwave 
field. The contribution to resistivity from trajectories in 
the bulk is analyzed in the frame of scattering on many 
well separated disk impurities. Here again the synchro- 
nization of cyclotron phase with the microwave phase 
takes place approximately at the same resonant val- 
ues. At these jV values the synchronization leads to long 
time capture of trajectories in disk vicinity. During this 
long time an initial cyclotron phase is washed out by 
noisy fluctuations and many rotations around disk and 
thus an electron escapes from a disk with an average zero 
displacement along the applied dc-field even if dynamics 
in disk vicinity is dissipative. This provides the main 
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mechanism of suppression of dissipative resistivity con- 
tribution from isolated impurities in the bulk. As a result 
the contribution of bulk to dissipative conductivity a^x 
is suppressed, as it was assumed in [3], and the main 
contribution to current is given by electron propagation 
along edge states stabilized by a microwave field. 

As we showed above the resonance width or resonance 
energy scale Ej. are approximately the same for the disk 
and wall cases (see Eqs. (|4j [S])). We note that for the 
disk case the energy Er is not sensitive to the disk radius 
as soon as it is significantly smaller than the cyclotron 
radius. Thus we expect that at jV- values the conductivity 
axx in the bulk is suppressed by a microwave field and 
at these fields the current is flowing essentially along sta- 
bilized edge states. In the case of Corbino geometry we 
have radial conductivity Urr which is determined by the 
bulk scattering and now the minima of (!„ are located at 



J. 



values (see e.g Figs. 12 13|14 where R 



CX (Tt 



The ZRS experiments performed in the Corbino geome- 
try give minima of arr at these jr values (see e.g. [^[?7] 1 
being in agreement with the synchronization theory. 

It is interesting to note that the nonlinear dynamics 
in vicinity of edge and disk impurity is well described by 
the Chirikov standard map |13| . The map description 
explains the location of resonances at integer values of j 
with an additional shift 5j « 1/4 produced by a finite 
separatrix width of nonlinear resonance. A finite rigidity 
of wall or disk potential can give a modification of this 
shift 5j. 

Our results show that the ZRS phases at appear 
only at weak noise corresponding to high mobility sam- 
ples. Strong noise destroys synchronization and trajecto- 
ries are no more captured at edge or disk vicinity. We also 
note that internal sample potentials with significant gra- 
dients act like a strong local rfc-field which destroys sta- 
bility regions around disk impurities or near edge. Due to 
that the ZRS effect exists only in high mobility samples. 
The resistance at ZRS minima drops significantly with 
the growth of microwave field strength since it increases 
the amplitude of nonlinear resonance which captures the 
synchronized trajectories. 

The synchronization theory of ZRS is based on classi- 
cal dynamics of noninteracting electrons. It is possible 
that electron-electron interaction effects can also sup- 
press the contribution to resistivity from neutral short 
range range scatterers (interface roughness, adatoms,...). 
Indeed, long capture times can increase the electron den- 
sity around these short ranged impurities transforming 
them into long range charged scatters that the other elec- 



trons can circumvent by adiabatically following the long 
range component of the disorder potential thereby avoid- 
ing a scattering event. However, the theoretical descrip- 
tion of this short-ranged impurity cloaking mechanism 
for the ZRS effect remains a serious challenge. 

Another important step remains the development of 
a quantum synchronization theory for ZRS. Even if in 
experiments the Landau level is relatively high v ^ 60, 
there are only about ten Landau levels inside a nonlin- 
ear resonance |S] and quantum effects should play a sig- 
nificant role. The general theoretical studies show that 
the phenomenon of quantum synchronization persists at 
small effective values of Planck constant Tief f but it be- 
comes destroyed by quantum fiuctuations at certain large 
values of h^f / [28 . 

The importance of quantum ZRS theory is also re- 
lated to a short range nature of the impurities consid- 
ered here, typically on a scale of a few nanometers or 
even less. We have modeled these impurities by disks 
with a radius that was only several times (in fact j 
times) smaller than the cyclotron radius which is not so 
close to microscopic reality. We could argue that in the 
quantum case a nanometer sized impurity would act ef- 
fectively as an impurity of a size of quantum magnetic 
length as ~ J'c/ ~ ?'c/8 ~ lOOnm. This gives a ratio 
fc/O'B ^ 8 which is comparable with the one used in our 
simulations with r^/r^ = j 7 but of course a quantum 
treatment of scattering on nanometer size impurity in 
crossed electric and magnetic and also microwave fields 
remains a theoretical challenge. We note that such type 
of scattering can be efficiently analyzed by tools of quan- 
tum chaotic scattering (see e.g. [211 [30]) and we expect 
that these tools will allow to make a progress in the quan- 
tum theory development of striking ZRS phenomenon. 

We hope that the synchronization theory of microwave 
induced ZRS phenomenon described here can be tested 
in further ZRS experiments. 
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